function res=And(lambda_t,x_t,n,xi,chi,mu_l,alpha,v,mu_x,rho,phi,psi,omega,eta,ngridpoints)

b=fliplr(psi*(1-omega.^(0:n-1))/(1-omega)-eta);

res=compute_exp(lambda_t,x_t,n,b,xi,chi,mu_l,alpha,v,mu_x,rho,phi,ngridpoints);